Experimental generation of steering odd dark beams of finite length 
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In this work we report the first realization of odd dark beams of finite length under controllable 
initial conditions. The mixed edge-screw phase dislocations are obtained by reproducing binary 
computer-generated holograms. Two effective ways to control the steering of the beams are analyzed 
experimentally and compared with numerical simulations. 
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I. INTRODUCTION 

Physically, Dark Spatial Solitons (DSSs) are localized 
intensity dips existing on stable background beams as 
a result of an exact counterbalance of diffraction and 
nonlinearity. A necessary condition for their existence 
is the presence of a phase dislocation in the wavefront 
along which the phase is indeterminate and the field am- 
plitude is zero. Besides the intriguing physical picture, 
the particular interest in the DSSs is motivated by their 
ability to induce gradient optical waveguides in bulk self- 
defocusing nonlinear media |l],g|3|,[|||. The only known 
truly two-dimensional (2D) DSSs are the Optical Vortex 
Solitons (OVSs) whereas in one transverse spatial di- 
mension DSSs manifest themselves as dark stripes ||. 
The odd initial conditions required to generate a funda- 
mental ID DSS correspond to an abrupt 7r-phase jump 
centered along the irradiance minimum of the stripe. The 
OVSs have a more complicated phase profile described by 
exp(imip) , where tp is the azimuthal coordinate in a plane 
perpendicular to the background beam propagation di- 
rection and m - the so-called topological charge (TC) - 
is an integer number. This phase function ensures a n- 
phase jump in each diametrical slice through the vortex 
core. Fundamental DSSs of these types have the common 
feature of zero transverse velocity if no perturbations are 
present. In contrast to that, ring dark solitary waves 
JiJ slowly change their parameters even when born from 
ideal odd initial conditions ||. 

In the pioneering work of Nye and Berry j^] it is con- 
jectured that mixed edge-screw dislocations cannot ex- 
ist. Despite that, almost two decades later an indica- 
tion for their existence was found |L0| for two interacting 
optical vortices of opposite topological charges. In our 
recent experiments on the generation of quasi-2D DSSs 
|nj we found that moderate saturation of the nonlinear- 
ity can stabilize the snake instability which usually leads 



to their decay. This made possible the first identification 
of ID Odd Dark Beams (ODBs) of finite length with their 
characteristic edge-screw phase dislocations [ pd]Jr2] | . The 
mixed dislocation forces the dark beams to steer in space. 
This appears to be of practical interest jl^] provided that 
there are effective ways to control the ODBs transverse 
velocity. 

In this article we report the first experimental real- 
ization of steering odd dark beams of finite length with 
mixed phase dislocations under controllable initial condi- 
tions. Two approaches to control their transverse veloc- 
ity are investigated experimentally and compared with 
numerical simulations. 



II. EXPERIMENTAL SETUP AND RESULTS 

A. Computer-generated holograms 

The phase portrait of the mixed edge-screw dislocation 
(see Fig. 1 in Ref. consists of a pair of semi- helices 
with a phase difference of ir to which an effective topo- 
logical charge of ±1/2 can be ascribed. Their spatial 
offset 2b determines the length of the edge part of the 
dislocation and ensures a phase jump of Aip in the di- 
rection perpendicular to the dark stripe of finite length. 
The phase function of this mixed phase dislocation can 
be described by 



§a,p{x, y) = A<pl arctan — — — 

7T \y + bf3 



(1-a) 



sgn(x) 
(2.1) 



where x and y denote the transverse Cartesian coordi- 
nates perpendicular and parallel to the dark beam. 2b 
stands for the length of the edge part of the dislocation 
and 
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f for \y\ < b 

a = I 1 and = -1 for y > b . (2.2) 
I 1 and (3 = 1 for y < — 

The pattern of the Computer-Generated Holograms 
(CGHs) used to produce this phase distribution consists 
of parallel lines which become curved at the position of 
the semi-vortex cores. In the edge part of the disloca- 
tion they terminate and reappear shifted, for a 7r-jump 
by one half of the pattern period. Holograms with such 
structures correspond to interference lines shifted along 
an imaginary line of finite length and to curved lines 
limiting the dislocations as observed in our previous ex- 
periment (see Fig. 8 in The binary CGHs used 
are photolithographically fabricated with a grating pe- 
riod of 20/im. Several holograms with various lengths 
of the edge part of the dislocation are etched on a com- 
mon substrate. Special attention was paid to align the 
edge parts of the different dislocations correctly on their 
common substrate. The simplicity of varying the dislo- 
cation length and magnitude of the phase jump is the 
main advantage |l3| ] of the approach we have chosen. 
The diffraction efficiency at first orders is measured to 
be 9%, close to the theoretical 10% limit for binary holo- 
grams. The unavoidable quantization inaccuracy of 7r/24 
|14f for holograms of this type is negligible for measure- 
ments with phase jumps of Aip = 37r/4, n, and 57r/4 
which are presented in this work. 

B. Experimental setup 

The setup used is similar to that in our previous ex- 
periments (see Fig. 1 in Ref. pl| ). Briefly, the beam 
of a single-line Ar + laser (A = 488nm) is used to re- 
construct the CGHs. The first-diffraction-order beam 
with the phase dislocation nested in is filtered through 
a slit and is focused on the entrance of the 10cm long 
Nonlinear Medium (NLM). After the desired propaga- 
tion path length, the beam is deflected by a prism im- 
mersed in the nonlinear liquid and is projected directly 
on a Charge-Coupled Device (CCD) camera array with 
a resolution of 13/im. The NLM is ethylene glycol dyed 
with DODCI (Lambdachrome) to reach an absorption co- 
efficient of 0.107cm -1 . In a calibration measurement we 
generated ID DSSs by using CGHs of the type described 
in Sec. |. The soliton constant la 2 (i.e. the product 
of the background-beam intensity / and the square of 
the dark beam width a measured at the 1/e- level) was 
found to reach its asymptotically constant value for in- 
put powers of P^f « 33mW. It is known that thermally 
self-defocusing liquids are both nonlocal and saturable. 
Since the saturation of the nonlinearity is able to modify 
the ODBs transverse velocity and profile, we needed to 
estimate it and to account for it in our numerical calcu- 
lations. In an independent measurement we realized a 
self-bending scheme similar to that used in [jl5],[l6| . The 
asymmetry required was introduced by an intentional tilt 



of the prism immersed in the NLM, which resulted in dif- 
ferent nonlinear propagation path-lengths for the differ- 
ent parts of the background beam. The strength of the 
self-bending effect was measured in the near field. For an 
absorptive nonlocal medium the choice of a suitable satu- 
ration model is not trivial ( p^j , see also Sec. IV of ifLsf). 
We found a good fit for the experimental data with the 
equation Ay ~ 1/(1 + J// SQt ) 7 . Using it we estimated 
P sa t ~ lOOrnW^ and 7 = 3. In addition to the careful 
alignment of the CGHs on the substrate, the holograms 
were reproduced to achieve vertical dark beam steering, 
which is not sensitive to possible undesired weak hori- 
zontal self-deflection of the background beam. Changing 
the ODB parameters (length-to-width ratio and magni- 
tude of the phase jump) was performed by a strict hori- 
zontal translation of the substrate. The accuracy of the 
alignment was tested by checking for equal steering of the 
ODBs reproduced from two identical holograms placed at 
opposite ends of the series of aligned CGHs. In this work 
the ODBs are identified by the corresponding lengths 2b 
of the edge portions of the dislocations in CGH pixels 
(lpix. = 5/jto) as encoded in the holograms. The de- 
flection Ax of the dark beams is measured in units of 
CCD camera pixels. We estimate that measurements 
with an encoded dislocation length of b/h[im CGH pixels 
corresponds to a dislocation length-to-ODB width ratio 
(2b/ a) = l/10&/5/im in the numerical simulations (e.g. 
14 pix. dislocation length corresponds to b/a = 1.4 in 
the simulations). 

C. ODB steering vs. dislocation length 

All data presented in this subsection refer to it phase 
jumps across the edge parts of the mixed dislocations. In 
Fig. la we plot the deflection Ax of ODBs with different 
dislocation lengths for input powers of 1.1 mW (circles), 
33mW (triangles) and 61mW (squares). The data ob- 
tained at P = 1.1 mW refer to a linear regime of propa- 
gation. The results at 33mW and 61mW are extracted 
from the experimental pictures shown in Fig. lb and 
Fig. lc, respectively, which are recorded for a nonlinear 
propagation path length of z = 8.5cm. The general ten- 
dency of linear increase of the deflection with decreasing 
the dislocation length is clearly expressed (Fig. la, solid 
line). The shortest mixed phase dislocation encoded was 
only lpix. long. In this case the strong deviation from 
the linear dependence is caused by the annihilation of 
the semi-charges due to the shortening of the edge part 
of the dislocation. It will be shown later that this short- 
ening accelerates for higher input powers (intensities). 
For that reason, even the ODB with an initially lQpix. 
long phase dislocation appears gradually less deflected 
at P = QlmW as compared to the case of P = 33mW 
(Fig. la). In Fig. lb,c the thick solid lines are intended 
to denote the positions of the ODBs at the entrance of 
the NLM. Because ODB steering is present also in the 



2 



linear regime of propagation (Fig. la, dots) these posi- 
tions (with respect to the dark beam intensity minimum 
and trailing peak maximum) are identified by numeri- 
cal simulations for b/a = 2.5. The identification corre- 
sponds to an encoded dislocation length of 25pix., i.e. 
to the most right frame shown in Fig. lb,c. Somewhat 
surprising is the weak sensitivity of the ODB deflection 
vs. background beam power (intensity) far from TC an- 
nihilation. It can be intuitively understood by recalling 
the known interaction scenario of well-separated OVSs of 
opposite TCs |ll||2(| . In this case the attraction between 
the OVSs is negligible as compared to their translation 
as a pair. At a constant input power of 33mW we mea- 
sured the ODB deflection vs. nonlinear propagation path 
length (Fig. 2). As expected, the ODB with lOpix. long 
dislocation has higher transverse velocity than that one 
with 22pix. long dislocation. The linearity in the depen- 
dencies is also well pronounced. 




z, cm 

FIG. 2. Deflection Aa; vs. nonlinear propagation path 
length z for encoded dislocation lengths of Wpix. (circles) 
and 22pix. (squares). The dashed and the solid lines are the 
respective linear fits. (P = 33mW; Aip — it). 



D. Phase control of the ODB steering 
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FIG. 1. Deflection of odd dark beams of finite length vs. 
dislocation length at input powers of 1.7m W (dots), 33mW 
(triangles), and 67mW (squares) for z = 8.5cm and Aip — ir 
(a). The points are extracted from the frames shown in Fig. lb 
(33mIF) and Fig.lc (67mW). From left to right the frames 
correspond to b — 1,10,14,18,22, and 25pix., respectively. 
The thick solid line indicates the calculated position of the 
ODB at the entrance of the NLM for 6 = 25pix.. 



As a second possible way to control the dark beam de- 
flection, we considered the variations in the magnitude 
of the phase jump Aip across the edge part of the mixed 
dislocation. In Fig. 3 we compare the experimental de- 
pendencies Ax(Aip) \b=22pix. (squares) and Ax(b) \Aip=n 
(dots). The straight lines are the respective linear fits. 
Because a problem in encoding a larger set of phase 
jumps in the CGHs was recognized too late, we mea- 
sured the deflection Aa; at Aip = 37r/4, it, and 5ir/4 only. 
In view of that the linear fit of the phase dependence 
in Fig. 3 appears to be assailable. Its linearity, how- 
ever, is confirmed by numerical simulations (see Sec. III). 
The dependence of Ax on Ax(Aip) and Ax(b) has been 
plotted in one figure in order to underline the fact that 
it appears to be easier to deflect the ODB by control- 
ling the phase than by controlling the dislocation length. 
The measurements are performed at a constant power of 
33mIT. 

Dislocation length, pix. 

8 10 12 14 16 18 20 22 24 26 
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FIG. 3. ODB deflection vs. magnitude of the phase jump 
Aip (hollow squares) and the dislocation length as encoded 
in the CGHs (dots). The straight lines are linear fits. 
(P = 33mlT). 
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E. Power /Intensity dependencies 



The ability of the dark spatial solitons |2|||[|j5| (and 
the dark spatial waves [^l"l , [22"|| ) to induce gradient all- 
optical waveguides in bulk sclf-defocusing NLM origi- 
nates in the negative nonlinear correction to the linear 
refractive indices of the media. In view of that, the in- 
tensity dependencies remain undoubtedly of interest, de- 
spite of the low sensitivity of the ODB steering on the 
input power (intensity). In Fig. 4 we present experimen- 
tal data on the power dependence of the length of the 
edge portion of the mixed dislocation for A<p — tt and 
for two different lengths of \Apix. and 22pix. encoded in 
the holograms. It is easy to understand that the (mixed) 
phase dislocations do not remain sharp and of an un- 
changed magnitude when the (odd) dark beam steers 
p3| . The dislocation lengths are estimated by evaluating 
the respective longitudinal ODB slices at 5% of the back- 
ground beam intensity (i.e. at the actual noise level in 
the recorded frames). Generally, the dislocation length 
decreases monotonically with increasing the input power. 
Asymptotically, the dislocation flattens and disappears, 
provided that the ratio b/a is less than 2. As mentioned 
in fl2]| , at b/a ~ 4 the ODBs should be expected to bend 
due to the snake instability [[TtJ . In fact we observed such 
a behavior for ODBs with encoded dislocation lengths 
ranging from 4 to 6. A vortex-beam creation is recog- 
nized by the convergence of two neighboring interference 
lines in one. However, the vortices formed by this insta- 
bility remained with highly overlapping cores JT^p^ ] . 




o-l . 1 . 1 . 1 . 1 . 1 
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FIG. 4. Edge dislocation length vs. input background 
beam power for two different dislocation lengths encoded in 
the CGHs. (Atp = iv;z = 8.5cm) 
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FIG. 5. ODB width (a) and length (b) (at the l/e level) vs. 
input background-beam power for two different dislocation 
lengths encoded in the CGHs. (z = 8.5cm) 

In Fig. 5a, b we plot the measured ODB widths (a) 
and their lengths (b) at the 1 /e-level as a function of the 
background beam power. The initial lengths of the mixed 
phase dislocations are denoted in pixels as encoded on the 
respective CGHs, whereas the widths and lengths at the 
exit of the NLM are measured in units of CCD camera 
pixels. The strong decrease in both the ODB width and 
length up to 17mW (approximately P^ /2) is followed 
by an approximate recovering of the width and length at 
approximately P}® ■ At higher input powers both trans- 
verse quantities decrease but the tendency is slower as 
compared to the situation below Pj® /2 and stabilizes 
asymptotically at high saturation. As it will be discussed 
in the next section, the minima in the curves plotted 
in Fig. 5a,b result from the reshaping of the OD beams 
at the particular l/e intensity-level chosen for evalua- 
tion. The estimation there shows that the first rapid de- 
crease in both transverse dimensions of the ODBs does 
not correspond to a 'soliton constant' formation. Gener- 
ally speaking, the ODBs analyzed are no solitary waves 
in the widely adopted sense, since they do not survive, for 
instance, a collision with a second ODB steering in the 
opposite direction |l^]. Nevertheless, the narrowing in 
both transverse directions for higher powers (intensities) 
should improve their guiding ability when signal beams 
or pulses are to be transmitted inside the ODBs and de- 
flected in space. This feature will be addressed elsewhere. 
It is interesting to note that, independent of the length 
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of the mixed dislocations (26), the widths (a) of the dis- 
locations should be initially equal (in the near field be- 
hind the CGHs) , but appeared different near the entrance 
of the NLM. The estimation has shown a CGH-to-NLM 
distance of approximately 4 Raylcigh diffraction lengths 
with respect to the initial ODB width. The well pro- 
nounced separation between the curves in Fig. 5a should 
be attributed to the different two-dimensional diffraction 
at different initial ODB length-to-width ratios. 



III. NUMERICAL SIMULATIONS 

In our numerical simulations we tried to model the 
experimentally obtained dependencies by accounting for 
the estimated moderate saturation of the nonlin earit y 
(P}$ = ssmw, P sat = lOOmW at 7 = 3; see Sec. |n§. 
The (2 + I)-dimensional nonlinear evolution of the steer- 
ing ODB of finite length in the bulk homogeneous and 
isotropic NLM is described by the generalized nonlinear 
Schrodinger equation 



dE 



drf 



E 



L 



Diff 



\E\ 



L NL (1 + S |£|2)7 



E = 0, 
(3.1) 



where the transverse spatial coordinates are normalized 
to the initial dark-beam width (( = x/a,rj = y/a), and 
the propagation path length is expressed in Rayleigh 
diffraction lengths LDiff = ka 2 . Further, Ljvl = 
{kri2lo)~ 1 is the nonlinear length, k is the wavenumbcr 
inside the NLM, and Iq is the background beam inten- 
sity. The adopted correction for the nonlinear refractive 
index is 



An = n 2 \E\ 2 /(l + s\E\ 2 )\ 



(3.2) 



with s = I P sa t = 0.3. As it was done in |JH|, the 
slowly varying electric-field amplitude of the ODB was 
chosen tan/i-shaped 



E(x,y) = V^oB(ri i0 (a;,?/))tanh 



O i*c>,/3 0,2/) 



(3.3) 



where r a ,p{x,y) = \/x 2 + a(y + /36) 2 is the effective 
Cartesian/radial coordinate, ^ a ,/3 is the phase dist ribu - 
tion of the mixed edge-screw dislo cation (see Eq. 2.1), 
and a and (3 are given by Eq. 2.2. The width w of the 



super-Gaussian background beam 



B(r) = exp 




(3.4) 



is chosen to exceed at least 15 times the ODB lengths. 
The model equation (3.1) was solved numerically by the 
beam propagation method on a 1024 x 1024 grid. It 



should be mentioned, that the initial width a(z — 0) of 
the ODB of finite length was chosen to correspond to that 



of an infinite ID ODSS (a 



l sol 



const. I, 



-1/2, 



It was 



proven numerically that the ODB deflection is insensitive 
to the particular value of a, and Figs. 6-8 are generated 
under this assumption. Actually, the nonlinearity causes 
an appreciable reshaping of the beams, in particular in 
the first evolution stage when the ODB starts steering, 
see Fig. 2 in Ref. |l^]). In order to improve the similar- 
ity between the experimental data (Figs. 4 and 5) and 
the numerical results (Figs. 9 and 10), an initial ODB 
width twice as large as in the experiment is assumed. In 
Fig. 6 we plot the ODB deflection vs. b/a \ z= o for differ- 
ent input powers (intensities). All data presented in this 
section refer to a normalized propagation path length of 
2 = 4:Lnl, which corresponds to that estimated for the 
experiment. Nevertheless all calculations are carried out 
up to 10-Ljvl whereby no qualitative deviation from the 
tendencies discussed is seen. The linearity in the ODB 
deflection vs. b/a is well obeyed except for b/a > 2.2. 
The longer ODBs steer slower, bend slightly, and decay 
into pairs of vortex beams for b/a > 4. In the linear 
regime of propagation, the ODBs also deflect but the 
deflection is stronger at higher input powers/intensities. 
This is more pronounced for shorter dislocations. In 
the experimental data (Fig. la), this behavior is much 
weaker, rather the deflection remains within the experi- 
mental accuracy. Looking for an adequate explanation, 
in a series of simulations we checked that a ±30% inaccu- 
racy in estimating P sat results only in ±5% deviation in 
the ODB deflection at z = ALnl- The observed tendency 
of a decrease of the ODB steering velocities at increased 
saturation is well understood [^5| but seems insufficient 
in quantity. We attribute the absence of a well expressed 
power dependence in Fig. la to the NLM nonlocality. In 
a separate experiment it was estimated, that the nonlo- 
cality in this medium is negligible on a spatial scale of 
several hundred of micrometers only [JL8|. 
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FIG. 6. Calculated ODB steering vs. b/a for different input 
powers with A(p — n. 
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FIG. 7. ODB steering along the NLM for b/a = 1.0 and 2.2. 
Crosses - some 10% of the numerical data. Vertical dashed 
line - propagation distance z/Lnl = 4 corresponding to the 
experimental conditions at 33m W . 

In Fig. 7 we plot the ODB deflection vs. the nonlin- 
ear propagation path length z/Lnl for b/a = 1.0 and 
2.2. As in the previous figure the magnitude of the edge 
part of the phase dislocation is set to Aip = it. Ones 
the ODB starts steering, its transverse velocity remains 
constant (see Fig. 2). The longer ODBs with longer edge 
dislocations, however, emit dispersive waves in their first 
evolution stage (z < ILml)- This causes a 'delay' in the 
steering along the NLM (Fig. 7, lower curve). 



2.5 




0.5 ■ 1 ■ 1 ■ 1 ■ 1 ■ r- 1 

0.5 0.6 0.7 0.8 0.9 1.0 

Phase jump, 71-units 

FIG. 8. Deflection of ODBs vs. phase jump Aip for 
b/a = 2.2 and for different input powers. 

In Fig. 8 we present results obtained for the phase- 
dependent control of the ODB deflection at different in- 
put powers (intensities). In qualitative agreement with 
the experimental observation, the linear increase in the 
ODB steering with decreasing the magnitude of the phase 
jump down to Aip — 0.5tt is evident. The compari- 
son of Figs. 6 and 8 confirms the conclusion that at a 
fixed nonlinear propagation distance the phase-controlled 
ODB deflection is more efficient as compared to that by 
varying the b/a ratio. 
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FIG. 9. ODB width vs. input power for Aip — ir and 
b/a = 1.4. (Dashed curve: full width at the 1/e- intensity 
level; solid: FWHM). 

Fig. 9 is intended to clarify the origin of the non- 
monotonic power dependencies of the ODB widths (and 
lengths) observed at input powers below Pj^ (see Fig. 5). 
The solid line represents the ODB full width at half max- 
imum (FWHM), the dashed one the full width at the 
1/e-intensity level. In these simulations Aip = tt and 
b/a =1.4 are assumed. Qualitatively, we obtained the 
same curves also for b/a = 2.2 by accounting for the ini- 
tial free-space propagation in the experiment (from the 
CGH to the entrance of the NLM; z w 4L Diff ). The 
minimum in the ODB width evaluated at the 1/e level 
originates from the reshaping of the beam profile which is 
caused by the moderate saturation. A similar reshaping 
is reported in Q (see Figs. 2-4 therein). At b/a = 2.2 
the data obtained for the ODB length vs. input power 
(intensity) were found to be even more sensitive to the 
intensity level of evaluation. Because of the transverse 
steering of the ODBs of finite length the edge portions 
of the dislocations shorten monotonically with increasing 
the background beam power (intensity) (Fig. 10). The 
numerical results are in a very good qualitative agree- 
ment with the experimental ones (Fig. 4). 




Power, mW 

FIG. 10. Length of the edge portion of the mixed phase 
dislocation vs. input power for Aip | z= o= tt and b/a = 1.4 
and 2.2, respectively. 



6 



IV. CONCLUSION 

The results presented show that the inherent steering 
dynamics of odd dark beams of finite length can be effec- 
tively controlled by varying both the magnitude Atp and 
the relative length b/a of the mixed edge-screw phase 
dislocation. The background-beam intensity has weak 
influence on the steering but is important for keeping 
the optically-induced gradient waveguides steep, which 
is crucial for all-optical guiding, deflection and switching 
of signal beams or pulses. Since the mixed phase dis- 
locations shorten and flatten along the nonlinear media 
(tending asymptotically to washout) the ODBs seem to 
be promising primarily for future short-range all-optical 
switching devices. 

ACKNOWLEDGMENTS 

A. D. would like to thank the Alexander von Humboldt 
Foundation for the award of a fellowship and the oppor- 
tunity to work in the stimulating atmosphere of the Max- 
Planck-Institut fur Quantenoptik (Garching, Germany). 
This work was also supported by the National Science 
Foundation of Bulgaria. 



[13] N. R. Heckenberg, R. McDuff, C. P. Smith, A. G. White, 

Opt. Lett. 17, 221 (1992). 
[14] W. -H. Lee, Prog. Optics XVI, 119 (1978). 
[15] A. E. Kaplan, JETP Lett. 9, 33 (1969). 
[16] M. S. Bordin, A. M. Kamuz, JETP Lett. 9, 351 (1969). 
[17] V. Tikhonenko, J. Christou, B. Luther-Davies, Yu. S. 

Kivshar, Opt. Lett. 21, 1129 (1996). 
[18] A. Dreischuh, G. G. Paulus, F. Zacher, F. Grasbon, H. 

Walther, Phys. Rev. E60, 6111 (1999). 
[19] K. Staminas, Chaos, Solitons & Fractals 4, 1783 (1994). 
[20] D. Neshev, A. Dreischuh, M. Assa, S. Dinev, Opt. Com- 

mun. 151, 413 (1998). 
[21] A. Dreischuh, V. Kamenov, S. Dinev, Appl. Phys. B 63, 

145 (1996). 

[22] D. Neshev, A. Dreischuh, S. Dinev, L. Windholz, J. Opt. 

Soc. Am. B 14, 2869 (1997). 
[23] R. Thurston, A. Weiner, J. Opt. Soc. Am. B 8, 471 (1991). 
[24] A.V. Mamaev, M. Saffman, A. A. Zozulya, Phys. Rev. Lett. 

76, 2262 (1996). 
[25] W. Krolikowski, N. Akhmediev, B. Luther-Davies, Phys. 

Rev. E 48, 3980 (1993). 
[26] V. Fikhonenko, Yu. S. Kivshar, V. Steblina, A. Zozulya, 

J. Opt. Soc. Am. B 15, 79 (1998). 



[1] Yu. S. Kivshar, B. Luther-Davies, Physics Reports 289, 81 

(1997) . 

[2] G. Swartzlander, Jr., C. Law, Phys. Rev. Lett. 69, 2503 
(1992). 

[3] E. A. Ostrovskaya, Yu. S. Kivshar, Opt. Lett. 23, 1268 

(1998) . 

[4] A.H. Carlsson, J.N. Malmberg, E.A. Ostrovskaya, T.J. 
Alexander, D. Anderson, M. Lisak, and Yu. Kivshar, lanl, 



patt-sol/9909009t 30 Sept. (1999). 
[5] C.T. Law, X. Zhang, and GA. Swartzlander, Opt. Lett. 
25, 55 (2000). 

[6] G. Allan, S. Skiner, D. Andersen, A. Smirl, Opt. Lett. 16, 
156 (1991). 

[7] Yu. S. Kivshar, X. Yang, Phys. Rev. E 50, R40-R43 (1994); 

ibid. Chaos, Solitons & Fractals 4, 1745 (1994). 
[8] D. Neshev, A. Dreischuh, V. Kamenov, I. Stefanov, S. 

Dinev, W. Fliefier, L. Windholz, Appl. Phys. B 64, 429 

(1997). 

[9] J. F. Nye, M. V. Berry, Proc. R. Soc. London A 336, 165 
(1974). 

[10] V. Bazhenov, M. Soskin, M. Vasnetsov, J. Mod. Opt. 39, 
985 (1992); I. Basistiy, V. Bazhenov, M. Soskin, M. Vas- 
netsov, Opt. Commun. 103, 422 (1993). 

[11] A. Dreischuh, G. G. Paulus, F. Zacher, Appl. Phys. B 69, 
113 (1999). 

[12] A. Dreischuh, G. G. Paulus, F. Zacher, I. Velchev, Appl. 
Phys. B 69, 107 (1999). 



7 



